My host's enemy is my enemy: plasmids carrying CRISPR-Cas as a defence against phages.

Bacteria are infected by mobile genetic elements like plasmids and virulent phages, and those infections significantly impact bacterial ecology and evolution. Recent discoveries reveal that some plasmids carry anti-phage immune systems like CRISPR-Cas, suggesting that plasmids may participate in the coevolutionary arms race between virulent phages and bacteria. Intuitively, this seems reasonable as virulent phages kill the plasmid's obligate host. However, the efficiency of CRISPR-Cas systems carried by plasmids can be expected to be lower than those carried by the chromosome due to continuous segregation loss, creating susceptible cells for phage amplification. To evaluate the anti-phage protection efficiency of CRISPR-Cas on plasmids, we develop a stochastic model describing the dynamics of a virulent phage infection against which a conjugative plasmid defends using CRISPR-Cas. We show that CRISPR-Cas on plasmids provides robust protection, except in limited parameter sets. In these cases, high segregation loss favours phage outbreaks by generating a population of defenceless cells on which the phage can evolve and escape CRISPR-Cas immunity. We show that the phage's ability to exploit segregation loss depends strongly on the evolvability of both CRISPR-Cas and the phage itself.

Forecasting antimicrobial resistance evolution.

Antimicrobial resistance (AMR) is a major global health issue. Current measures for tackling it comprise mainly the prudent use of drugs, the development of new drugs, and rapid diagnostics. Relatively little attention has been given to forecasting the evolution of resistance. Here, we argue that forecasting has the potential to be a great asset in our arsenal of measures to tackle AMR. We argue that, if successfully implemented, forecasting resistance will help to resolve the antibiotic crisis in three ways: it will (i) guide a more sustainable use (and therefore lifespan) of antibiotics and incentivize investment in drug development, (ii) reduce the spread of AMR genes and pathogenic microbes in the environment and between patients, and (iii) allow more efficient treatment of persistent infections, reducing the continued evolution of resistance. We identify two important challenges that need to be addressed for the successful establishment of forecasting: (i) the development of bespoke technology that allows stakeholders to empirically assess the risks of resistance evolving during the process of drug development and therapeutic/preventive use, and (ii) the transformative shift in mindset from the current praxis of mostly addressing the problem of antibiotic resistance a posteriori to a concept of a priori estimating, and acting on, the risks of resistance.

estimateR: an R package to estimate and monitor the effective reproductive number.

BACKGROUND: Accurate estimation of the effective reproductive number ([Formula: see text]) of epidemic outbreaks is of central relevance to public health policy and decision making. We present estimateR, an R package for the estimation of the reproductive number through time from delayed observations of infection events. Such delayed observations include confirmed cases, hospitalizations or deaths. The package implements the methodology of Huisman et al. but modularizes the [Formula: see text] estimation procedure to allow easy implementation of new alternatives to the currently available methods. Users can tailor their analyses according to their particular use case by choosing among implemented options. RESULTS: The estimateR R package allows users to estimate the effective reproductive number of an epidemic outbreak based on observed cases, hospitalization, death or any other type of event documenting past infections, in a fast and timely fashion. We validated the implementation with a simulation study: estimateR yielded estimates comparable to alternative publicly available methods while being around two orders of magnitude faster. We then applied estimateR to empirical case-confirmation incidence data for COVID-19 in nine countries and for dengue fever in Brazil; in parallel, estimateR is already being applied (i) to SARS-CoV-2 measurements in wastewater data and (ii) to study influenza transmission based on wastewater and clinical data in other studies. In summary, this R package provides a fast and flexible implementation to estimate the effective reproductive number for various diseases and datasets. CONCLUSIONS: The estimateR R package is a modular and extendable tool designed for outbreak surveillance and retrospective outbreak investigation. It extends the method developed for COVID-19 by Huisman et al. and makes it available for a variety of pathogens, outbreak scenarios, and observation types. Estimates obtained with estimateR can be interpreted directly or used to inform more complex epidemic models (e.g. for forecasting) on the value of [Formula: see text].

The effect of combining antibiotics on resistance: A systematic review and meta-analysis.

When and under which conditions antibiotic combination therapy decelerates rather than accelerates resistance evolution is not well understood. We examined the effect of combining antibiotics on within-patient resistance development across various bacterial pathogens and antibiotics. We searched CENTRAL, EMBASE and PubMed for (quasi)-randomised controlled trials (RCTs) published from database inception to November 24th, 2022. Trials comparing antibiotic treatments with different numbers of antibiotics were included. A patient was considered to have acquired resistance if, at the follow-up culture, a resistant bacterium was detected that had not been present in the baseline culture. We combined results using a random effects model and performed meta-regression and stratified analyses. The trials' risk of bias was assessed with the Cochrane tool. 42 trials were eligible and 29, including 5054 patients, were qualified for statistical analysis. In most trials, resistance development was not the primary outcome and studies lacked power. The combined odds ratio (OR) for the acquisition of resistance comparing the group with the higher number of antibiotics with the comparison group was 1.23 (95% CI 0.68-2.25), with substantial between-study heterogeneity (I2 =77%). We identified tentative evidence for potential beneficial or detrimental effects of antibiotic combination therapy for specific pathogens or medical conditions. The evidence for combining a higher number of antibiotics compared to fewer from RCTs is scarce and overall, is compatible with both benefit or harm. Trials powered to detect differences in resistance development or well-designed observational studies are required to clarify the impact of combination therapy on resistance.

Sex and gender in infection and immunity: addressing the bottlenecks from basic science to public health and clinical applications.

Although sex and gender are recognized as major determinants of health and immunity, their role is rarely considered in clinical practice and public health. We identified six bottlenecks preventing the inclusion of sex and gender considerations from basic science to clinical practice, precision medicine and public health policies. (i) A terminology-related bottleneck, linked to the definitions of sex and gender themselves, and the lack of consensus on how to evaluate gender. (ii) A data-related bottleneck, due to gaps in sex-disaggregated data, data on trans/non-binary people and gender identity. (iii) A translational bottleneck, limited by animal models and the underrepresentation of gender minorities in biomedical studies. (iv) A statistical bottleneck, with inappropriate statistical analyses and results interpretation. (v) An ethical bottleneck posed by the underrepresentation of pregnant people and gender minorities in clinical studies. (vi) A structural bottleneck, as systemic bias and discriminations affect not only academic research but also decision makers. We specify guidelines for researchers, scientific journals, funding agencies and academic institutions to address these bottlenecks. Following such guidelines will support the development of more efficient and equitable care strategies for all.

A Genome-Scale Metabolic Model of Marine Heterotroph Vibrio splendidus Strain 1A01.

While Vibrio splendidus is best known as an opportunistic pathogen in oysters, Vibrio splendidus strain 1A01 was first identified as an early colonizer of synthetic chitin particles incubated in seawater. To gain a better understanding of its metabolism, a genome-scale metabolic model (GSMM) of V. splendidus 1A01 was reconstructed. GSMMs enable us to simulate all metabolic reactions in a bacterial cell using flux balance analysis. A draft model was built using an automated pipeline from BioCyc. Manual curation was then performed based on experimental data, in part by gap-filling metabolic pathways and tailoring the model's biomass reaction to V. splendidus 1A01. The challenges of building a metabolic model for a marine microorganism like V. splendidus 1A01 are described. IMPORTANCE A genome-scale metabolic model of V. splendidus 1A01 was reconstructed in this work. We offer solutions to the technical problems associated with model reconstruction for a marine bacterial strain like V. splendidus 1A01, which arise largely from the high salt concentration found in both seawater and culture media that simulate seawater.

Estimation and worldwide monitoring of the effective reproductive number of SARS-CoV-2.

The effective reproductive number Re is a key indicator of the growth of an epidemic. Since the start of the SARS-CoV-2 pandemic, many methods and online dashboards have sprung up to monitor this number through time. However, these methods are not always thoroughly tested, correctly placed in time, or are overly confident during high incidence periods. Here, we present a method for timely estimation of Re, applied to COVID-19 epidemic data from 170 countries. We thoroughly evaluate the method on simulated data, and present an intuitive web interface for interactive data exploration. We show that, in early 2020, in the majority of countries the estimated Re dropped below 1 only after the introduction of major non-pharmaceutical interventions. For Europe the implementation of non-pharmaceutical interventions was broadly associated with reductions in the estimated Re. Globally though, relaxing non-pharmaceutical interventions had more varied effects on subsequent Re estimates. Our framework is useful to inform governments and the general public on the status of epidemics in their country, and is used as the official source of Re estimates for SARS-CoV-2 in Switzerland. It further allows detailed comparison between countries and in relation to covariates such as implemented public health policies, mobility, behaviour, or weather data.

Over the past two and a half years, countries around the globe have struggled to control the transmission of the SARS-CoV-2 virus within their borders. To manage the situation, it is important to have an accurate picture of how fast the virus is spreading. This can be achieved by calculating the effective reproductive number (Re), which describes how many people, on average, someone with COVID-19 is likely to infect. If the Re is greater than one, the virus is infecting increasingly more people, but if it is smaller than one, the number of cases is declining. Scientists use various strategies to estimate the Re, which each have their own strengths and weaknesses. One of the main difficulties is that infections are typically recorded only when people test positive for COVID-19, are hospitalized with the virus, or die. This means that the data provides a delayed representation of when infections are happening. Furthermore, changes in these records occur later than measures that change the infection dynamics. As a result, researchers need to take these delays into account when estimating Re. Here, Huisman, Scire et al. have developed a new method for estimating the Re based on available data records, statistically taking into account the above-mentioned delays. An online dashboard with daily updates was then created so that policy makers and the population could monitor the values over time. For over two years, Huisman, Scire et al. have been applying their tool and dashboard to COVID-19 data from 170 countries. They found that public health interventions, such as mask requirements and lockdowns, did help reduce the Re in Europe. But the effects were not uniform across the globe, likely because of variations in how restrictions were implemented and followed during the pandemic. In early 2020, the Re only dropped below one after countries put lockdowns or other severe measures in place. The Re values added to the dashboard over the last two years have been used pro-actively to inform public health policies in Switzerland and to monitor the spread of SARS-CoV-2 in South Africa. The team has also recently released programming software based on this method that can be used to track future disease outbreaks, and extended the method to estimate the Re using SARS-CoV-2 levels in wastewater.

The effect of sequencing and assembly on the inference of horizontal gene transfer on chromosomal and plasmid phylogenies.

The spread of antibiotic resistance genes on plasmids is a threat to human and animal health. Phylogenies of bacteria and their plasmids contain clues regarding the frequency of plasmid transfer events, as well as the co-evolution of plasmids and their hosts. However, whole genome sequencing data from diverse ecological or clinical bacterial samples are rarely used to study plasmid phylogenies and resistance gene transfer. This is partially due to the difficulty of extracting plasmids from short-read sequencing data. Here, we use both short- and long-read sequencing data of 24 clinical extended-spectrum [Formula: see text]-lactamase (ESBL)-producing Escherichia coli to estimate chromosomal and plasmid phylogenies. We compare the impact of different sequencing and assembly methodologies on these phylogenies and on the inference of horizontal gene transfer. We find that chromosomal phylogenies can be estimated robustly with all methods, whereas plasmid phylogenies have more variable topology and branch lengths across the methods used. Specifically, hybrid methods that use long reads to resolve short-read assemblies (HybridSPAdes and Unicycler) perform better than those that started from long reads during assembly graph generation (Canu). By contrast, the inference of plasmid and antibiotic resistance gene transfer using a parsimony-based criterion is mostly robust to the choice of sequencing and assembly method. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Turnover in Life-Strategies Recapitulates Marine Microbial Succession Colonizing Model Particles.

Particulate organic matter (POM) in the ocean sustains diverse communities of bacteria that mediate the remineralization of organic complex matter. However, the variability of these particles and of the environmental conditions surrounding them present a challenge to the study of the ecological processes shaping particle-associated communities and their function. In this work, we utilize data from experiments in which coastal water communities are grown on synthetic particles to ask which are the most important ecological drivers of their assembly and associated traits. Combining 16S rRNA amplicon sequencing with shotgun metagenomics, together with an analysis of the full genomes of a subset of isolated strains, we were able to identify two-to-three distinct community classes, corresponding to early vs. late colonizers. We show that these classes are shaped by environmental selection (early colonizers) and facilitation (late colonizers) and find distinctive traits associated with each class. While early colonizers have a larger proportion of genes related to the uptake of nutrients, motility, and environmental sensing with few pathways enriched for metabolism, late colonizers devote a higher proportion of genes for metabolism, comprising a wide array of different pathways including the metabolism of carbohydrates, amino acids, and xenobiotics. Analysis of selected pathways suggests the existence of a trophic-chain topology connecting both classes for nitrogen metabolism, potential exchange of branched chain amino acids for late colonizers, and differences in bacterial doubling times throughout the succession. The interpretation of these traits suggests a distinction between early and late colonizers analogous to other classifications found in the literature, and we discuss connections with the classical distinction between r- and K-strategists.

Epidemiological and evolutionary consequences of different types of CRISPR-Cas systems.

Bacteria have adaptive immunity against viruses (phages) in the form of CRISPR-Cas immune systems. Currently, 6 types of CRISPR-Cas systems are known and the molecular study of three of these has revealed important molecular differences. It is unknown if and how these molecular differences change the outcome of phage infection and the evolutionary pressure the CRISPR-Cas systems faces. To determine the importance of these molecular differences, we model a phage outbreak entering a population defending exclusively with a type I/II or a type III CRISPR-Cas system. We show that for type III CRISPR-Cas systems, rapid phage extinction is driven by the probability to acquire at least one resistance spacer. However, for type I/II CRISPR-Cas systems, rapid phage extinction is characterized by an a threshold-like behaviour: any acquisition probability below this threshold leads to phage survival whereas any acquisition probability above it, results in phage extinction. We also show that in the absence of autoimmunity, high acquisition rates evolve. However, when CRISPR-Cas systems are prone to autoimmunity, intermediate levels of acquisition are optimal during a phage outbreak. As we predict an optimal probability of spacer acquisition 2 factors of magnitude above the one that has been measured, we discuss the origin of such a discrepancy. Finally, we show that in a biologically relevant parameter range, a type III CRISPR-Cas system can outcompete a type I/II CRISPR-Cas system with a slightly higher probability of acquisition.

Applying mixture model methods to SARS-CoV-2 serosurvey data from Geneva.

Serosurveys are an important tool to estimate the true extent of the current SARS-CoV-2 pandemic. So far, most serosurvey data have been analyzed with cutoff-based methods, which dichotomize individual measurements into sero-positives or negatives based on a predefined cutoff. However, mixture model methods can gain additional information from the same serosurvey data. Such methods refrain from dichotomizing individual values and instead use the full distribution of the serological measurements from pre-pandemic and COVID-19 controls to estimate the cumulative incidence. This study presents an application of mixture model methods to SARS-CoV-2 serosurvey data from the SEROCoV-POP study from April and May 2020 in Geneva (2766 individuals). Besides estimating the total cumulative incidence in these data (8.1% (95% CI: 6.8%-9.9%)), we applied extended mixture model methods to estimate an indirect indicator of disease severity, which is the fraction of cases with a distribution of antibody levels similar to hospitalized COVID-19 patients. This fraction is 51.2% (95% CI: 15.2%-79.5%) across the full serosurvey, but differs between three age classes: 21.4% (95% CI: 0%-59.6%) for individuals between 5 and 40 years old, 60.2% (95% CI: 21.5%-100%) for individuals between 41 and 65 years old and 100% (95% CI: 20.1%-100%) for individuals between 66 and 90 years old. Additionally, we find a mismatch between the inferred negative distribution of the serosurvey and the validation data of pre-pandemic controls. Overall, this study illustrates that mixture model methods can provide additional insights from serosurvey data.

Estimating plasmid conjugation rates: A new computational tool and a critical comparison of methods.

Plasmids are important vectors for the spread of genes among diverse populations of bacteria. However, there is no standard method to determine the rate at which they spread horizontally via conjugation. Here, we compare commonly used methods on simulated and experimental data, and show that the resulting conjugation rate estimates often depend strongly on the time of measurement, the initial population densities, or the initial ratio of donor to recipient populations. Differences in growth rate, e.g. induced by sub-lethal antibiotic concentrations or temperature, can also significantly bias conjugation rate estimates. We derive a new 'end-point' measure to estimate conjugation rates, which extends the well-known Simonsen method to include the effects of differences in population growth and conjugation rates from donors and transconjugants. We further derive analytical expressions for the parameter range in which these approximations remain valid. We present an easy to use R package and web interface which implement both new and previously existing methods to estimate conjugation rates. The result is a set of tools and guidelines for accurate and comparable measurement of plasmid conjugation rates.

Constrained optimization of divisional load in hierarchically organized tissues during homeostasis.

It has been hypothesized that the structure of tissues and the hierarchy of differentiation from stem cell to terminally differentiated cell play a significant role in reducing the incidence of cancer in that tissue. One specific mechanism by which this risk can be reduced is by minimizing the number of divisions-and hence the mutational risk-that cells accumulate as they divide to maintain tissue homeostasis. Here, we investigate a mathematical model of cell division in a hierarchical tissue, calculating and minimizing the divisional load while constraining parameters such that homeostasis is maintained. We show that the minimal divisional load is achieved by binary division trees with progenitor cells incapable of self-renewal. Contrary to the protection hypothesis, we find that an increased stem cell turnover can lead to lower divisional load. Furthermore, we find that the optimal tissue structure depends on the time horizon of the duration of homeostasis, with faster stem cell division favoured in short-lived organisms and more progenitor compartments favoured in longer-lived organisms.

Test-trace-isolate-quarantine (TTIQ) intervention strategies after symptomatic COVID-19 case identification.

The test-trace-isolate-quarantine (TTIQ) strategy, where confirmed-positive pathogen carriers are isolated from the community and their recent close contacts are identified and pre-emptively quarantined, is used to break chains of transmission during a disease outbreak. The protocol is frequently followed after an individual presents with disease symptoms, at which point they will be tested for the pathogen. This TTIQ strategy, along with hygiene and social distancing measures, make up the non-pharmaceutical interventions that are utilised to suppress the ongoing COVID-19 pandemic. Here we develop a tractable mathematical model of disease transmission and the TTIQ intervention to quantify how the probability of detecting and isolating a case following symptom onset, the fraction of contacts that are identified and quarantined, and the delays inherent to these processes impact epidemic growth. In the model, the timing of disease transmission and symptom onset, as well as the frequency of asymptomatic cases, is based on empirical distributions of SARS-CoV-2 infection dynamics, while the isolation of confirmed cases and quarantine of their contacts is implemented by truncating their respective infectious periods. We find that a successful TTIQ strategy requires intensive testing: the majority of transmission is prevented by isolating symptomatic individuals and doing so in a short amount of time. Despite the lesser impact, additional contact tracing and quarantine increases the parameter space in which an epidemic is controllable and is necessary to control epidemics with a high reproductive number. TTIQ could remain an important intervention for the foreseeable future of the COVID-19 pandemic due to slow vaccine rollout and highly-transmissible variants with the potential for vaccine escape. Our results can be used to assess how TTIQ can be improved and optimised, and the methodology represents an improvement over previous quantification methods that is applicable to future epidemic scenarios.

Ecological effects of stress drive bacterial evolvability under sub-inhibitory antibiotic treatments.

Stress is thought to increase mutation rate and thus to accelerate evolution. In the context of antibiotic resistance, sub-inhibitory treatments could then lead to enhanced evolvability, thereby fuelling the adaptation of pathogens. Combining wet-lab experiments, stochastic simulations and a meta-analysis of the literature, we found that the increase in mutation rates triggered by antibiotic treatments is often cancelled out by reduced population size, resulting in no overall increase in genetic diversity. A careful analysis of the effect of ecological factors on genetic diversity showed that the potential for regrowth during recovery phase after treatment plays a crucial role in evolvability, being the main factor associated with increased genetic diversity in experimental data.

Plasmid co-infection: linking biological mechanisms to ecological and evolutionary dynamics.

As infectious agents of bacteria and vehicles of horizontal gene transfer, plasmids play a key role in bacterial ecology and evolution. Plasmid dynamics are shaped not only by plasmid-host interactions but also by ecological interactions between plasmid variants. These interactions are complex: plasmids can co-infect the same cell and the consequences for the co-resident plasmid can be either beneficial or detrimental. Many of the biological processes that govern plasmid co-infection-from systems that exclude infection by other plasmids to interactions in the regulation of plasmid copy number-are well characterized at a mechanistic level. Modelling plays a central role in translating such mechanistic insights into predictions about plasmid dynamics and the impact of these dynamics on bacterial evolution. Theoretical work in evolutionary epidemiology has shown that formulating models of co-infection is not trivial, as some modelling choices can introduce unintended ecological assumptions. Here, we review how the biological processes that govern co-infection can be represented in a mathematical model, discuss potential modelling pitfalls, and analyse this model to provide general insights into how co-infection impacts ecological and evolutionary outcomes. In particular, we demonstrate how beneficial and detrimental effects of co-infection give rise to frequency-dependent selection on plasmid variants. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Increase in antimicrobial resistance in Escherichia coli in food animals between 1980 and 2018 assessed using genomes from public databases.

BACKGROUND: Next-generation sequencing has considerably increased the number of genomes available in the public domain. However, efforts to use these genomes for surveillance of antimicrobial resistance have thus far been limited and geographically heterogeneous. We inferred global resistance trends in Escherichia coli in food animals using genomes from public databases. METHODS: We retrieved 7632 E. coli genomes from public databases (NCBI, PATRIC and EnteroBase) and screened for antimicrobial resistance genes (ARGs) using ResFinder. Selection bias towards resistance, virulence or specific strains was accounted for by screening BioProject descriptions. Temporal trends for MDR, resistance to antimicrobial classes and ARG prevalence were inferred using generalized linear models for all genomes, including those not subjected to selection bias. RESULTS: MDR increased by 1.6 times between 1980 and 2018, as genomes carried, on average, ARGs conferring resistance to 2.65 antimicrobials in swine, 2.22 in poultry and 1.58 in bovines. Highest resistance levels were observed for tetracyclines (42.2%-69.1%), penicillins (19.4%-47.5%) and streptomycin (28.6%-56.6%). Resistance trends were consistent after accounting for selection bias, although lower mean absolute resistance estimates were associated with genomes not subjected to selection bias (difference of 3.16%±3.58% across years, hosts and antimicrobial classes). We observed an increase in extended-spectrum cephalosporin ARG blaCMY-2 and a progressive substitution of tetB by tetA. Estimates of resistance prevalence inferred from genomes in the public domain were in good agreement with reports from systematic phenotypic surveillance. CONCLUSIONS: Our analysis illustrates the potential of using the growing volume of genomes in public databases to track AMR trends globally.

Pathogen invasion-dependent tissue reservoirs and plasmid-encoded antibiotic degradation boost plasmid spread in the gut.

Many plasmids encode antibiotic resistance genes. Through conjugation, plasmids can be rapidly disseminated. Previous work identified gut luminal donor/recipient blooms and tissue-lodged plasmid-bearing persister cells of the enteric pathogen Salmonella enterica serovar Typhimurium (S.Tm) that survive antibiotic therapy in host tissues, as factors promoting plasmid dissemination among Enterobacteriaceae. However, the buildup of tissue reservoirs and their contribution to plasmid spread await experimental demonstration. Here, we asked if re-seeding-plasmid acquisition-invasion cycles by S.Tm could serve to diversify tissue-lodged plasmid reservoirs, and thereby promote plasmid spread. Starting with intraperitoneal mouse infections, we demonstrate that S.Tm cells re-seeding the gut lumen initiate clonal expansion. Extended spectrum beta-lactamase (ESBL) plasmid-encoded gut luminal antibiotic degradation by donors can foster recipient survival under beta-lactam antibiotic treatment, enhancing transconjugant formation upon re-seeding. S.Tm transconjugants can subsequently re-enter host tissues introducing the new plasmid into the tissue-lodged reservoir. Population dynamics analyses pinpoint recipient migration into the gut lumen as rate-limiting for plasmid transfer dynamics in our model. Priority effects may be a limiting factor for reservoir formation in host tissues. Overall, our proof-of-principle data indicates that luminal antibiotic degradation and shuttling between the gut lumen and tissue-resident reservoirs can promote the accumulation and spread of plasmids within a host over time.

Correction: Practical considerations for measuring the effective reproductive number, Rt.

[This corrects the article DOI: 10.1371/journal.pcbi.1008409.].

Multidrug Resistance Dynamics in Salmonella in Food Animals in the United States: An Analysis of Genomes from Public Databases.

The number of bacterial genomes deposited each year in public databases is growing exponentially. However, efforts to use these genomes to track trends in antimicrobial resistance (AMR) have been limited thus far. We used 22,102 genomes from public databases to track AMR trends in nontyphoidal Salmonella in food animals in the United States. In 2018, genomes deposited in public databases carried genes conferring resistance, on average, to 2.08 antimicrobial classes in poultry, 1.74 in bovines, and 1.28 in swine. This represents a decline in AMR of over 70% compared to the levels in 2000 in bovines and swine, and an increase of 13% for poultry. Trends in resistance inferred from genomic data showed good agreement with U.S. phenotypic surveillance data (weighted mean absolute difference ± standard deviation, 5.86% ± 8.11%). In 2018, resistance to 3rd-generation cephalosporins in bovines, swine, and poultry decreased to 9.97% on average, whereas in quinolones and 4th-generation cephalosporins, resistance increased to 12.53% and 3.87%, respectively. This was concomitant with a decrease of blaCMY-2 but an increase in blaCTX-M-65 and gyrA D87Y (encoding a change of D to Y at position 87). Core genome single-nucleotide polymorphism (SNP) phylogenies show that resistance to these antimicrobial classes was predominantly associated with Salmonella enterica serovar Infantis and, to a lesser extent, S. enterica serovar Typhimurium and its monophasic variant I 4,[5],12:i:-, whereas quinolone resistance was also associated with S. enterica serovar Dublin. Between 2000 and 2018, trends in serovar prevalence showed a composition shift where S. Typhimurium decreased while S. Infantis increased. Our findings illustrate the growing potential of using genomes in public databases to track AMR in regions where sequencing capacities are currently expanding. IMPORTANCE Next-generation sequencing has led to an exponential increase in the number of genomes deposited in public repositories. This growing volume of information presents opportunities to track the prevalence of genes conferring antimicrobial resistance (AMR), a growing threat to the health of humans and animals. Using 22,102 public genomes, we estimated that the prevalence of multidrug resistance (MDR) in the United States decreased in nontyphoidal Salmonella isolates recovered from bovines and swine between 2000 and 2018, whereas it increased in poultry. These trends are consistent with those detected by national surveillance systems that monitor resistance using phenotypic testing. However, using genomes, we identified that genes conferring resistance to critically important antimicrobials were associated with specific MDR serovars that could be the focus for future interventions. Our analysis illustrates the growing potential of public repositories to monitor AMR trends and shows that similar efforts could soon be carried out in other regions where genomic surveillance is increasing.